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SUMMARY 

A method of analyzing the first mode stability and unbalance response of 
multimass flexible rotors is presented whereby the multimass system is modeled 
as an equivalent single mass modal model including the effects of rotor flex- 
ibility, general linearized hydrodynamic journal bearings, squeeze film bearing 
supports and rotor aerodynamic cross coupling. Expressions for optimum bearing 
and support damping are presented for both stability and unbalance response. 

The method is Intended to be used as a preliminary design tool to quickly 
ascertain the effects of bearing and support changes on rotor-bearing system 
performance. 


INTRODUCTION 

The purpose of this paper is to present a method that is easily applied 
to a large class of industrial turbomachines which will predict the instabil- 
ity onset speeds including the effects of shaft flexibility, generalized hydro- 
dynamic journal bearings, squeeze film bearing supports and aerodynamic rotor 
excitation. By restricting the analysis to the class of turbomachines which 
have their mass centers inboard of the bearings and are relatively symmetric 
about the mass center, a simplified rotor model is established that allows the 
equations of motion to be manipulated analytically. The restriction to turbo- 
machines of this class does not severely limit its usefulness or applicability 
since the majority of industrial machines fall within this class. A major 
advantage of the method is that it allows machine designers to quickly deter- 
mine the effects of bearing changes, shaft modifications and bearing support 
designs to determine appropriate system designs. Those designs deemed appro- 
prxate using the method can be further verified and finalized using more gen- 
eral analysis tools. A large time and cost savings should be realized by 
using the method to eliminate impractical designs without incurring the high 
computer costs and large amounts of data reduction time required using more 
general rotor dynamic analysis programs. The method is intended to supplement 
the more general techniques, not to replace them. 
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Specifically the method represents the rotor, bearings and supports as a 
linear dynamical system with the rotor represented by a generalized modal mass 
and modal stiffness. The Inclusion of aerodynamic effects using the Alford 
model shows that an optimum bearing or support damping exists which will maxi- 
mize the aerodynamic forces required to produce rotor instability for a given 
rotor and bearing or supfort flexibility. Definite expressions for the optunum 
damping are derived which eliminates the requirement of a parametric variation 
to find the optimum damping, 

A number of instability mechanisms have been identified over the years. 

Chief among these are hydrodynamic bearings and aerodynamic effects. 

(1) Hydrodynamic journal bearings - Most industrial turbomachine rotors 
are supported in fluid film bearings. Bearings with fixed geometries develop 
forces within the lubricant film which couple the motion of the rotor in any ^ 
two orthogonal transverse directions and can produce instability under certain 
conditions of rotor speed and bearing load. Recent emphasis has been placed on 
determining the linearized fluid film force coefficients for various bearing 
geometries which can be used for linear stability analyses (refs. 1-9). 

(2) Interaction with working fluid flow (Aerodynamic cross coupling) 

The vibratory motion of a rotor within the machine working fluid can produce 
differential pressures around the rotor which result in destabilizing forces 
and moments. Because of the very complicated nature of the flows around tur- 
bine blades and through centrifugal impeller passages a general analysis 
describing the working fluid forces as a function of rotor motion has not been 
developed. The coupling between the working fluid and the rotor is commonly 
referred to as aerodynamic cross coupling. 

The most commonly used approximation for predicting the magnitude of aero- 
dynamic cross coupling assumes that the radial deflection of the rotor produces 
a variation -'n the thermodynamic efficiency around the impeller circumference 
and hence a torque variation (ref. 10). The differential torque produces a 
force on the impeller proportional to the displacement and perpendicular to ^ 
in the direction of impeller rotation and thus tends to drive the rotor in 
a forward (same direction as rotation) whirl. The cross coupling force on each 
impeller is a function of stage mean torque, pitch radius of the impeller, blade 
length and impeller displacement. A damping force is not postulated.^ The force 
is also proportional to a generally unknown constant which is a function of 
fluid mass flow rate, pressure, enthalpy and other fluid properties (ref. 10). 

Figure 1 shows the frequency spectra for a 7-stage industrial centrifugal 
compressor designed to operate at 13500 rpm. Originally mounted in rigidly 
supported tilting pad bearings, the machine became unstable when the rotor 
speed exceeded 10500 rpm. The instability was due to aerodynamic cross coupling 
and seal forces. The spectra in figure la illustrates the subsynchronous fre- 
quency representing the instability shown by the 4300 cpm component of the 
spectra. This frequency was identified as the first damped frequency of the 
system and remained nearly constant as the machine speed was increased. The 
synchronous unbalance excitation is indicated by the IN line and is seen to 
peak at about 4300 rpm. The vibration at the subsynchronous frequency was 
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ufficiently large to pervent sustained operation at 13500 rpm and the machine 
operation was restricted to below 10000 rpm for a considerable period of tLe 
Monetary losses to the user due to reduced production capability were large. 

Figure Ih shows the spectra after installation of properly designed 
Sstabn tilting pad bearing locations. The subsynchronous 

opelat^ elix,l„ated. The unit haa been aucLssLur 

on^sShUlS (ref U)“" capability for nearly four years with no indication 
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ANALYSIS 

Equations of Motion 

Consider the generic rotor-bearing system shown schematically in figure 2. 
This system is representative of a large class of rotors supported in two bear- 
ings V7ith most of the mass between the bearings. These machines generally op- 
erate above the first critical speed and the instability that occurs is gener- 
ally associated with the first mode of the system (ref. 12). The homogeneous 
equations of motion for the system are 


mx + k{x - x)-hqy *0 
2 s 2 X ^ 


( 1 ) 
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where q represents a destabilizing cross-coupled aerodynamic force acting at 
the rotor center. 

It is convenient to nondimensionalize the equations of motion using the 
following variables 
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where c is the bearing radial clearance and O) is the rotor speed. Substitution 
into equations (l)-(4) yields the following set of nondimensional equations 
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where the primes denote differentiation with respect to cot. 
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This set of linear ordinary differential equations may be solved to deter- 
mine the stability of the rotor-bearing system subject to the bearing forces 
and the equivalent aerodynamic force. The solution to these equations has the 
form 
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i = 1, 2 


where s represents the system complex eigenvalues, 
s = P + iA 

The stability of the system is determined by the sign of P. If P < 0, the 
system is stable and small motions about the system equilibrium position will 
decay with time. If P > 0, the system is pronounced unstable and the motions 
will Increase with time. If P = 0, the system is said to be at the onset of 
instability or instability threshold. It is this condition that will be exam- 
ined in greater detail. For any given set of conditions, A represents the 
rotor whirl ratio, that is, the ratio of the rotating shaft damped natural 
frequency to the rotor speed. 

For any given combination of linearized bearing characteristics and aero- 
dynamic excitation, there is a rotor speed where the real part of the eigen- 
value, P, is zero. This rotor speed is the instability onset speed, and the 
rotor is stable for speeds below this. 

Since P = 0, by definition, at the instability onset speed, the solution 
to the equations of motion becomes 
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where is the instability onset speed for the given bearing and aerodynamic 
conditions. Substitution into equations (4)-(8) gives 
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Since equations (9) and (10) are real, they may easily be solved for X 
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_ A nontrivial solution to this set of equations exists only if the deter- 
minant of the matrix of coefficients of X and Y is zero. Since the coeffi- 
cients are complex, the determinant is alio complex. The condition of a zero 
etermxnant requires that the real and imaginary parts of the determinant 
both be zero. Expanding equations (15) and (16) into matrix form and solving 
for the real and imaginary parts of the determinant results in the following 
two equations. 
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Equations (17) and (18) represent two equations in the two unknows (QX) ^ and 
. Once these variables are known, the Instability onset speed for the given 
conditions can be detericined and is 

n - ^ 

s X 


or in dimensional terms 
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( 22 ) 


The solution is obtained in the following way. First, write equation (18) 
as 
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For hydrodynamic journal bearings, - C (ref, 5) so the last term in 

equation (23) is zero and is a funttion^only of the known bearing coeffi- 
cients. 
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Now equation (20) is used to determine D^ and equation (17) is rearranged to 
give X^ . 
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Finally, equation (21) or (22) is used to determine the instability onset 
speed. The proper use of equation (25) requires some explanation. If the 
entire term within the radical is positive, equation (26) gives two real value 
for (QX)^ and, hence, two instability onset speeds are predicted. For the giv> 
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bearing and aerodynamic conditions, the rotor will be stable only for rotor 
speeds between the two predicted onset speed values. This condition occurs 
only for the case when the aerodynamic excitation is greater than zero. As 0 
approaches zero, the lower value of O) approaches zero and the rotor will be 
Stable for all rotor speeds below the upper value of co . 

s 

_ It is apparent from examination.of the radical in equation (25) that for 
a given level of shaft flexibility, 6, and a given set of bearing coefficients 
which uniquely determine D, , there is an upper bound to the level of aerodynamic 
excitation, Q, that will allow the terms inside the radical to be positive! 

Once this value of Q is surpassed, the radical term becomes imaginary. This 
implies that the rotor will be unstable at any rotor speed. Thus, Imposing the 
condition that the radical term be zero in equation (25) will yield the maximum 
value of aerodynamic excitation for which a stable operating speed exists. 

Therefore, the maximum aerodynamic excitation the rotor can withstand occurs 
when 
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Rearranging gives 
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Bearing Induced Instability 
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to examine in detail the effect of hydrodynamic bearings 
fl^^ible rotors. With aerodynamic excitation Q = 0, 
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Noting that with Q - 0, the whirl ratio at the instability onset speed is inde- 
pendent of shaft flexibility, 6, the rigid rotor onset speed is 
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Substitution into equation (30) gives the flexible rotor instability onset 
speed as 


n 

Q = (33) 

® /l + 

which clearly shows the reduction in stability with increasing shaft flexibility. 
For a given flexibility parameter, 6, the rotor instability onset speed can be 
increased by choosing bearings with a higher bearing instability threshold speed, 
n , and lower whirl ratio, X. As becomes very large, the limiting value 

of ^instability onset speed becomes 
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or in dimensional terms 
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Hence the instability onset speed of a flexible rotor is proportional to the 
rigid bearing natural frequency of the rotor. From equations (24) and (26) it 
is observed that increased bearing principal stiffness and damping and decreased 
cross-coupled stiffness increase bearing induced instability onset speeds. 
Further increases may be obtained from asymmetry in the principal stiffness and 
damping . 


Aerodynamically Induced Instability 

If aerodynamic cross-coupling is present it was previously noted that a 
maximum value of Q exists beyond which a particular flexible rotor-bearing 
system will be unstable at all speeds. That value was shown to be 


^ 6(1 + SDi) 


With Q = 0 , equation (25) becomes 
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Substituting into equation (20) gives the effective cross coupling at the 
bearings due to aerodynamic excitation 
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(38) 
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This value of is then used to calculate . 
The instability onset speed becomes 
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It should be remembered that equation (38) represents the instability onset 
speed with maximum aerodynamic cross-coupling present. With smaller values of 
Q the actual onset speed may be higher or lower as will be subsequently shown. 

Example 1 Single-Stage Centrifugal Pump in Plain Cylindrical Journal Bearings . 

A single-stage centrifugal pump has the impeller centrally mounted on the shaft. 
The rotor has the following characteristics: 

Shaft Length = 107,32 cm 
Shaft Diameter = 10.16 cm 
Impeller Weight = 1800,0 N 
Operating Speed = 4,000 rpm 

The rotor is mounted in two identical plain cylindrical journal bearings 
with the following characteristics: 

Bearing Length = 5.08 cm 
Bearing Diameter = 10.16 cm 
Radial Clearance = 0,0508 mm 
Oil Viscosity = 1.23 N-s/cm^ 

The shaft weight is 667,2 N. Hence, the modal weight is approx im.ately 
half the shaft weight plus the impeller weight, or 2134.5 N. The effective 
shaft stiffness is = 420,283 N/cm which corresponds to 6 = 1.0. The bearing 
coefficients are 
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The operating speed parameter is 
= 0-953 

0 O^g 

Equations (24), (36), (38), and (39) give the maximum aerodynamic excitation as 
0^ = 0.344 
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and the instability onset speed as 

a = 1.57 
s 

Since ^ ^ , the rotor will be unstable with Q present. The variation of 
instability onset speed with Q is shown by the solid curve in figure 3. This 
figure was obtained by varying Q is equation (25). The allowable aerodynamic 
excitation is Q = 0.17, The dashed curve in figure 3 shows the effect of 
changing the bearing length from 5.08 cm to 4.45 cm which increases the allow- 
able aerodynamic excitation from 0.17 to 0.245. 


Optimum Damping 

The stability analysis presented in the preceding sections show that for 
a given set of nondimensional bearing coefficients, the permissible aerodynamic 
rotor excitation is a function of rotor operating speed and a maximum value is 
indicated . 

The fact that a maximum value of allowable aerodynamic excitation exists 
may be construed as the existence of an optimum damping. This concept may be 
visualized in the following way. 


Since the nondimensional stiffness and damping coefficients are constant 
for a given operating condition, increasing the nondimensional speed (the 
ordinate of stability map) is equivalent to decreasing the dimensional damping 
in the bearings since damping is inversely proportional to speed. The stability 
map may be thought of as the effect of increasing the rotor speed assuming the 
nondimensional bearing stiffness and damping coefficients are independent of 
speed. The optimum dimensional damping occurs when and is 
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where C and C are the actual nondimensional damping coefficients for the 
operatifti condiHons being considered and is the dimensional speed at 

i.e. , 



This value of 9, is obviously chosen as the optimum since the maximum 
allowable aerodynamic excitation, Qjj^, occurs at this value. For values of 
Q < Q , the actual dimensional damping is excessive, and for 9 > the 

dampin|™is insufficient. The optimum nondimensional damping coefficients are, 
therefore, 
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If the bearings or supports are assumed symmetric, an explicit expression 
for optimum damping may be obtained. Letting 
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the following relationships are obtained 
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The instability onset speed becomes 
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For the system to tolerate the maximum aerod3m.amic excitation, the operating 

speed must be the same as ^ , 

sm* 


1 + 6(K^ + K^) 


V 


5(1 + 25k ) 

r 

2(1 + 


5\) 


(45) 


0 


= Q 

sm 


and the optimum damping is 
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(46) 
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It may be shown (ref . 13 , 14 , 15) that the effective damping acting at the 
rotor center is 


C = 
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where a is the ratio of actual bearing damping to optimum bearing damping, i.e., 


a = 
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and C is the effective damping with optimum bearing damping. Since the 
allowIBle aerodynamic is proportional to the effective damping, an estimate 
of the permissible aerodynamic excitation with non-optimum damping is 


Q = 



(48) 


Equation (46) can also be written in terms of the bearing (support) prin- 
cipal stiffness to shaft stiffness ratio 
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where K = 2 


Flexible, Damped Bearing Supports 

The equations previously presented may also be used to evaluate the effects 
of squeeze film bearing supports in series with a hydrodynamic bearing. Such 
a support is shown in figure 4 (ref. 16). In this particular application, the 
squeeze film is used to support a ball bearing mounted rotor, although any type 
of shaft bearing may be utilized. The particular values of the squeeze film 
bearing parameters are, of course, dependent on the rotor-bearing system it 
will support. The essential features of the squeeze film bearing include an 
annular clearance between the squeeze film journal and damper housing, lubri- 
cant supply grooves or holes, bearing end seals, and journal mechanical retainer 
springs. Various combinations of these components are used to achieve appro- 
priate stiffness and damping properties for a particular application. Since 
the squeeze film journal does not rotate, the hydrodynamic squeeze film forces 
result only from translational motion of the journal which squeezes the lubri- 
cant . 


Figure 5 schematically illustrates a flexible, damped bearing-support com- 
ponent. The bearing and suppport are represented by general linearized stiff- 
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ness and damping force coefficients. In many squeeze film bearing applications, 
the inner bearing mass, m^, is statically supported in mechanical springs to 
align and preload the inner bearing within the squeeze film annulus. Therefore, 
the general support radial stiffness is represented by hydrodynamic stiffness 
coefficients, mechanical retainer spring stiffness 

coefficients, bearing and support can be combined to give 

overall effective stiffness and damping properties (refs. 2, 13). However, 
care must be exercised when doing so. The effective coefficients are functions 
of the whirl ratio, X. Therefore, the use of the stability equations requires 
an iterative process whereby a value of X is chosen, the effective coefficients 
are calculated, and the value of X predicted by equation (26) is compared with 
the chosen value. The process is repeated until the calculated and assumed 
values of X are the same. A judicious first guess for X can be made since the 
damped natural frequency will usually lie somewhere between the undamped natural 
frequency of the rotor-bearing-support system and the undamped rigid bearing, 
rigid support natural frequency. Also, the whirl ratio dependence will be much 
less if either the bearings or supports have a large impedance due to high stiff- 
ness or damping. 

A further consideration is that under some conditions the effective damping 
matrix will not be symmetric (c ^ ^ c . With this condition, D , given by 

equation (23), is a function of^tli^ aerodynamic cross coupling of thh system as 
well as the bearing coefficients. Equations (20), (23), and (26) must be solved 

iteratively for D and D . 

1 2 

Example 2 11-Stage Centrifugal Compressor 

An 11— stage centrifugal compressor is shown in figure 6, The rotor physical 
characteristics are 

Overall Length = 216 cm 
Bearing Span = 175.3 cm 

Rotor Weight = 5782.4 N 

Operating Speed = 10,000 rpm 

The first rigid bearing critical speed is N = 3894 rpm. The second critical 
speed is above the 10,000 rpm operating speed with the bearings described above. 
The rotor is very nearly symmetrical about the bearing midspan and gyroscopic 
effects are negligible. The rotor modal weight is 2891.2 N, and the shaft stiff- 
ness is k = 490350 N/cm. 
s 

The rotor was originally mounted in tilting pad bearings with the following 
dimensions : 

Length = 5.72 cm 

Diameter = 11.43 cm 

Radial Pad Clearance = 0.1016 mm 

Bearing Radial Clearance = 0.0508 mm 

Number of Pads = 5 

Pad Arc Length = 60° 

Offset Factor = 0.50 
Preload = 0.50 
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Load Direction = Pad Pivot 

Oil Viscosity = 1.38 x 10“^ N-s/cm^ 

The bearing stiffness and damping coefficients are shown in figure 7 (ref. 7). 

At 10000 rpm the bearing coefficients are = K = 36.23, C^x ~ ^yy ~ ^3.63. 
Field experience indicated that the level of aerodynamic excitation at 10000 
rpm was q = 35000 N/cm and rotor instability was observed. 

Applying equations (36)-(39) the maximum permissible aerodynamic excitation 
Q = 0.078 (c] = 11098 N/cm). This is considerably less than the anticipated 
value of q = 35OT0 N/cm. Since occurs at $7^^ = 1.50 and the operating speed 

er is ~ ^ i even less aerodynamic excitation can occur . The vari~ 
of with Q is shown in figure 8. For ~ 3.37, the allowable 

aerodynamic excitation is Q = 0.055 (q = 7826 N^cm) . 

If the bearing clearance and pad clearance are Increased and pad clearance 
are increased to 0.1143 mm and 0.2286 mm respectively, the bearing coefficients 
at 10,000 rpm are 

K, = 8.83 K = 14.86 

= 9.88 = 13.29 

For this condition = 0.970 (q^ = 61340 N/cm) at ^ ~ 1.40. The operating 

speed parameter is = 5.06 (the change is due to the change in clearance. 

The variation of with Q for this case is also shown in figure 8. The allow- 
able value of aerodynamic excitation is Q = 0.59 (q = 37310 N/cm). This is 
slightly larger than the anticipated value. Therefore, it is apparent that 
large changes in stability can be achieved simply by modifying bearing clearance 

CONCLUSIONS 

(1) The single mass representation including generalized linearized bearing 
and aerodynamic excitation yields equations which predict instability onset 
speeds due to bearing cross coupling and aerodynamic effects. 

(2) Expressions for the optimum bearing damping for given bearing stiffness 
coefficients which will maximize the allowable aerodynamic excitation at a given 
operating speed and predict the maximum aerodynamic excitation have been devel- 
oped . 

(3) The form of the optimum damping expressions for stability are applic- 
able to first mode synchronous unbalance response and the optimum damping mini- 
mizes that response. 

(4) Optimum nondimensional bearing damping is a function of the rotr oper- 
ating speed, rigid bearing undamped critical speed, bearing clearance and ratio 
of principal bearing stiffness to the fundamental rotor modal. 

(5) The effect of flexible, damped bearing supports can be incorporated 
into the single mass modal representation using equivalent bearing-support 
linear stiffness and damping coefficients. 
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FREQUENCY SPECTRUM FilO’tCPMt 


FREQUENCY SPECTRUM Fr)O'(CPM) 


(a) (b) 

Pig. 1. Frequency Spectra of 7-Stngc Centrifugal Compresaor 

(a) Subaynctironous InsCnhlliCy with N > 10500 RPM 

(b) Stabilization with Vfcll Dcoigned Squeeze I’llm Bearing 



Fig. 2 


Schematic of Single Mass Representation 
of Symmetric Multimass Rotors 

















